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Direct numerical simulation (DNS) and linear analysis (LIA) of isotropic tur- 
bulence interacting with a shock wave are performed for several upstream shock 
normal Mach numbers (Mi ). Turbulence kinetic energy (TKE) is amplified across 
the shock wave, but this amplification tends to saturate beyond Mj = 3.0. TKE 
amplification and Reynolds stress anisotropy obtained in DNS are consistent with 
LIA predictions. Rapid evolution of TKE immediate downstream of the shock wave 
persists for all shock strengths and is attributed to the transfer between kinetic and 
potential modes of turbulence energy through acoustic fluctuations. Changes in en- 
ergy spectra and various length scales across the shock wave are predicted by LIA, 
which is consistent with DNS results. Most turbulence length scales decrease across 
the shock. Dissipation length scale (pq 1 * 3 /e), however, increases slightly for shock 
waves with Mi < 1.65. Fluctuations in thermodynamic variables behind the shock 
wave stay nearly isentropic for Mi < 1.2 and deviate significantly from isentropy 
for the stronger shock waves due to large entropy fluctuation generated through the 
interaction. 


1. Motivation and objective 

The presence of shock waves is an important feature that distinguishes high-speed 
supersonic flows. Understanding the mechanisms of turbulence interacting with a 
shock wave is not only of generic interest, but also of fundamental importance in 
predicting the interactions of turbulent boundary layers with the shock waves which 
occur in many engineering applications. Since the 1950’s, linear analyses (LIA) on 
the modification of elementary disturbance waves, such as vortical, acoustic, and 
entropic waves, by the shock wave have been performed with an emphasis on the 
acoustic wave generation behind the shock wave (Ribner 1953, 1954, 1968, Moore 
1953, Kerrebrock 1956, Chang 1957, McKenzie and Westphal 1968). Recently, the 
applicability of homogeneous Rapid Distortion Theory (RDT) on shock/turbulence 
interaction was investigated by Jacquin et al. (1993). 

There has been a significant accumulation of experimental data on the shock tur- 
bulence interaction during the last decade. Interaction of turbulent boundary layers 
with a shock wave over a corner was investigated by many research groups, among 
them are Dolling and Or (1985), Andreopoulos and Muck (1987), Smits and Muck 
(1987). A general finding from these experiments is that Reynolds shear stress and 
turbulence intensities are amplified across the shock wave. The studies of oblique 
shock wave/turbulent boundary layer interaction included several additional phe- 
nomena which complicated the flow behavior (Honkan and Andreopoulos 1992). To 
isolate the effects of a shock wave on turbulence, several experiments (Debieve and 
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Lachaxme 1986, Keller and Merzkirch 1990, Jacquin et al. 1991, Honkan et a l. 1992) 
on the interaction between the shock wave and grid-generated turbulence have been 
performed. They found that turbulence is amplified and turbulence length scales in- 
crease across the shock wave. But the length scale increase contradicts the intuitive 
expectation that mean flow compression should decrease the relevant turbulence 
length scales. The issue of length scale change will be thoroughly discussed in the 

present paper (Sec. 2.2). . . 

Numerical simulations of the shock turbulence interaction are just beginning to 
emerge. Using a shock capturing numerical technique, Rotman (1991) calculated 
the change in a two dimensional turbulent flow caused by the passage of the traveling 
shock wave. He found that the shock causes an increase in the turbulent kinetic 
energy and that the length scale of the turbulent field is reduced upon passage 
of the shock. Lee et al. (1991a, 1992) conducted direct numerical simulations of 
two and three dimensional turbulence interacting with a shock wave. They found 
that vorticity amplification compared well with the linear analysis predictions, and 
turbulent kinetic energy undergoes rapid increase behind the shock wave. The 
spectrum was found to be enhanced more at large wave numbers, leading to an 
overall length scale decrease. 

In the present report, interaction of isotropic turbulence with a strong shock wave 
is studied to investigate the effects of shock strength on turbulence modification. 
A numerical technique to simulate turbulence interacting with a strong shock wave 
without resolving its structure was developed, and it validated this technique against 
the shock-resolving simulations (Lee 1993). The simulation results are compared 
with the results from a linear analysis, and they are contrasted against the results 
from the weak shock case to show the shock strength effects. 


2. Accomplishments 

The parameters of the simulation are the mean Mach number (Mi), the fluctua- 
tion Mach number (M t ), and the turbulence Reynolds number based on the Taylor 
microscale ( Re x ) upstream of the shock wave. In the simulation, all of the turbu- 
lence scales are fully resolved, while the effect of the shock wave on turbulence is 
captured (rather than fully resolved). Two new simulations are conducted for the 
interaction with strong shock waves M\ = 2.0, 3.0), and the results from shock- 
resolving simulations (Lee et al. 1993) for the interaction with a weak shock wave 
— 1.05, 1.1, 1.2) tire quoted to investigate the effects of the shock normal Mach 
number. Table 1 lists the simulation parameters, where the values of M t and Re\ 
are taken at the location immediately upstream of the shock. 


Table 1. Parameters for the simulations of shock turbulence interaction 


Case 

Mi 

M t 

Re a 

Icq 

A 

2.0 ~~ i 

0.108 

19.0 

4.0 

B 

3.0 

0.110 

19.7 

4.0 
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FIGURE 1(a). Evolution of the normal components of the Reynolds stress: lines 

for M\ = 2.0, and symbols for M\ = 3.0. , • : i? u , , x ; R 22 , , 

+ : R33. Vertical lines denote the boundar ies of shock intermittency. 

2.1 Turbulence velocity fluctuation 

Interaction of turbulence with a shock wave generates acoustic waves downstream 
of the shock, part of which undergo rapid decay (Ribner 1953). LIA predicts that 
turbulent kinetic energy is amplified across the shock wave and the decaying acoustic 
waves contribute significantly to the streamwise fluctuations just behind the shock 
wave. 

Figure 1(A) shows the evolution of the diagonal components of the Reynolds stress 
tensor, R l} — u"u". The off-diagonal components stay close to zero over the entire 
flow field since turbulence is isotropic upstream and axisymmetric downstream of the 
shock. The streamwise component in the shock region contains the intermittency 
effects due to the oscillations of the shock. For more details of the intermittency 
effects on turbulence statistics, see Lee et al. (1992). The boundaries of the shock 
oscillations are defined as the locations where mean dilatation dui fdx t = 0; du\ jdx 1 
is negative inside the shock wave and slightly positive away from the shock due to 
viscous heating. All the velocity fluctuations are enhanced during the interaction. 
The velocity fluctuations are axisymmetric behind the shock, and their return to 
isotropy is negligible compared to the decay. Away from the shock wave, all the 
velocity fluctuations decay monotonically due to the viscous dissipation. 

Mach number dependence of the far-field velocity fluctuation amplification pre- 
dicted by LIA is shown in Figure 1(B). All components of the velocity fluctuation 
are amplified across the shock wave, and the amplification of TKE tends to satu- 
rate beyond Mi = 3.0. The shock normal component is amplified more for shock 
waves with Mi < 2.0 while the opposite is true for Mi > 2.0. In DNS, however, the 
streamwise velocity fluctuation away from the shock is larger than the transverse ve- 
locity fluctuations, which apparently contradicts with the LIA prediction. Viscous 




332 


Sang s an Lee 



FIGURE 1(b). Amplification of velocity fluctuation va rian ces acr oss the sho ck 
wave predicted by LIA at far away from the shock. u" 2 , u" 2 & u' 3 ' 2 . 



Figure 1(c). Evolution of velocity fluctuation variances behind the shock wave 
predicted by LIA (Mi = 2.0). u" 2 , u'J & uj . 

TKE dissipation rate behind the shock for the transverse components are found to 
be signifi can tly higher than the streamwise component. Therefore, comparing low 
Reynolds number DNS results directly with the inviscid linear analysis is not fair. 
After the viscous decay is compensated for by extrapolating the curves to the shock 
location, the trend of amplification from the DNS is found to be consistent with the 
LIA prediction. 
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The rapid evolution of velocity fluctuations which was observed for weak shock 
turbulence interaction (Lee et ai 1993) persists in the present simulations of strong 
shock turbulence interaction. In order to understand the downstream evolution 
of the velocity fluctuations, the budget of the Reynolds stress transport equation 
downstream of the shock wave is investigated. All terms in the transport equation 
can be accurately evaluated, since all the flow variables are fully resolved both 
in time and space outsi de the shock wave. As for weak shock case, the pressure 
transport term (-(p'u")^) in the inhomogeneous (or the shock normal) direction 
is mainly responsible for the rapid evolution of the streamwise velocity fluctuation. 
The evolution of the velocity fluctuations downstream of the shock wave predicted 
by LIA is shown in Figure 1(C), which reproduces the main feature of the rapid 
evolution from the DNS. Hence, the rapid TKE evolution behind the shock wave 
can be explained mainly as a linear process. This rapid evolution in the streamwise 
velocity fluctuation is due to a correlation between the vortical and decaying acoustic 
fluctuations behind the shock wave. The acoustic velocity fluctuations and vortical 
velocity fluctuations are anti-correlated just behind the shock, and the correlation 
between the two fluctuations decreases rapidly as the amplitude of the acoustic 
wave decays exponentially away from the shock wave. In previous studies (Lee et 
ai 1991a, 1992, 1993), the correlations between vortical and acoustic waves were 
not properly accounted for, and the prediction capability of the linear analysis was 
not fully appreciated. 

Another facet of the rapid evolution of velocity fluctuations is revealed by an 
equation for linear acoustic energy balance (Thompson 1985). The continuity and 
momentum equations for the linearized fluctuating components can be written as 


dp ' „ dp' 

dt Uk dxk 


+ P 


M 

dx * 


= 0, 


dul ^ dv£ 1 _ da" 

dt + k dx k +pdxi dx k 


by assuming that there exist no mean flow gradients, where pa" ( = r" ) denotes the 
viscous stress. For an ideal gas, an infinitesimal density fluctuation can be related 
to the pressure and entropy (s) fluctuations by 


fL _ _ fL 

p ip c? 

where c p is the specific heat at constant pressure. Multiplying the continuity equa- 
tion by p', contracting the momentum equations by u", and cancelling density- 
dilatation correlation by using above thermodynamic relation with neglecting en- 
tropy fluctuation effect (— s'u" -/c p c), the following equation (in the averaged form) 
follows. 


d 

dx k 


i p' u '1 

c ' 2c 2 2 p 2 7 pc 


c 3 dxk 


= 0 . 
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FIGURE 2(a). Evolution of the integrated quantities from the acoustic energy 

balance for M\ = 2.0: A, 5, C, A + B + C, A + B + 

CAD. 



k 0 x i 

FIGURE 2(b). Evolutions of pressure-dilatation correl ation a nd entropy-dilatation 

correlation for Mi =2.0: pV/ k hpu 0 k 0 , s’u'l k /c p u 0 k 0 . Vertical lines 

denote the boundaries of shock intermit tency. 


If this relation is satisfied the phenomenon can be explained in terms of linear 
acoustic energy balance. 

The acoustic balance equation is integrated in the streamwise direction from the 
downstream side of the shock (:r*) to give 
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U I u"u" 
c 2 c 2 


in 


+ 


Ui/ 
c 2 p 2 


i*i 


+ 


1 p'u" 

7 l>c 

V* 

C 


Xl /*• < da" , 

- / dr, -- 0, 

«. J*. c dx k 


where ||/||‘ = f(b) - f(a). The integrated results are shown in Figure 2(A). In all 
the cases we investigated (with different shock strengths and upstream turbulence 
intensities) the acoustic energy balance is satisfied with little deviation. The rapidly 
evolving acoustic energy — sum of scaled density and velocity fluctuations — is 
found to be mainly balanced by the pressure transport once the decay due to viscous 
dissipation is compensated. Therefore, the rapid evolution of velocity fluctuations 
can be attributed to the acoustic energy balance: energy transfer from the acoustic 
potential energy in the form of density (or pressure) fluctuations to turbulence 
kinetic energy. This is consistent with the fact that the pressure-transport term 
is scaled best by flow variables associated with acoustic wave propagation (Lee et 
al. 1993). Note that density fluctuation is replaced by the pressure fluctuation 
using the isentropic relation in deriving the acoustic energy balance equation, even 
though entropy fluctuation behind the shock wave contributes significantly to the 
density fluctuation (as is shown in Sec. 2.3). This is justified because as shown 
in Figure 2(B) the neglected entropy-dilatation correlation is found to be less than 
30% of the pressure-dilatation correlation in the zone of interest, and 5% of the 
pressure transport term. Even though thermodynamic fluctuations are far from 
isentropic, the contribution of entropy fluctuations to the acoustic energy balance 
can be neglected. The entropy-dilatation correlation vanishes in the linear limit 
and the acoustic energy balance derived above holds exactly in the linear analysis, 
which ignores viscous dissipation. 

Variance of vorticity fluctuation is a main contributor to the TKE dissipation 
rate. Figure 3(A) shows the evolution of vorticity components. The transverse 
components are amplified across the shock, while the streamwise component is 
hardly affected. Mach number dependence of transverse vorticity variance ampli- 
fication predicted by LIA is shown in Figure 3(B). LIA predicts no amplification 
of the streamwise component. The amplification trend and its amplification ratio 
obtained from DNS Eire found to be consistent with the LIA prediction. 


2.2 Turbulence length scales 

Experimental studies (Debieve et al. 1986, Keller et al. 1990, Honkan et al. 
1992) have reported that large scale turbulent motions are enhanced more than 
small scale motions as turbulence passes through a shock wave, leading to the 
overall increase of turbulence length scales, especially of microscales. LIA predicts 
that Taylor microscales decrease across the shock wave for all shock strengths, 
which was confirmed by DNS for weak shock waves (Lee et al. 1991a, 1993). For 
weak shock waves, chcinges in some turbulence length scEiles were too small to draw 
definite conclusions on the issue. 

To investigate the scEile-dependent amplification of turbulence, the modification 
of power spectra across the shock wave (M, = 2.0) is computed through LIA for 
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Figure 3(a). Evolution of_the vorticity fluctuation variances with Mi - 2.0: 



FIGURE 3(b) . Amplification of transverse vorticity fluctuation variances predicted 
by LIA. 

the one dimensional spectrum in the shock-normal (longitudinal) and transverse 
direction which is shown in Figure 4(A) and (B), respectively. In the longitudinal 
spectrum, significant scale- dependent amplification is observed: more amplification 
at small scales than at large scales. Large scale part of E 2 {h) is even suppressed 
through the interaction. In the transverse spectrum, more amplification at small 
scales is found for EM and E 2 (k 2 ), while more amplification at large scale is 
found for _E 3 (fc 2 ). The energy spectrum used in the analysis is the von karman 
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FIGURE 4. (a) LIA prediction of the shock-normal direction (vs. iti) one dim- 
ensional energy spectra change across the shock wave for M x = 2.0: E x : up- 
stream, downstream, E 2 & £3: upstream, downstream, (b) LIA 

prediction of the transverse direction (vs. k 2 ) one dimensional energy spectra change 

across the shock wave for M x = 2.0: E 2 : upstream, downstream, E x : 

upstream, downstream, £ 3 : upstream, downstream. 


spectrum (Hinze 1975), but the results obtained in the analysis are insensitive to 
choice of the spectrum. Since the spectrum amplification pattern is different for 
different spectrum (e.g. £i(A: 2 ) or £ 3 (Ar 2 )), the issue of the length scale change 
should be addressed for the specific length scale only. In the following, changes in 
various turbulence length scales are discussed. 

To directly check the scale-dependent turbulence amplification, transverse power 
spectra of velocity fluctuations in a numerically simulated field from case A are 
shown for upstream and downstream of the shock wave in Figure 5. Amplification is 
more pronounced at the large wave numbers, which is consistent with the prediction 
by the linear analysis in Figure 4(B). 

teller et al. (1990) reported that both the density microscale and the integral 
length scale in the shock normal direction increase for shock waves with M x < 1.24. 
In the present simulation, the spectrum changes of density and temperature fluc- 
tuations across the shock are found to be similar to those of velocity fluctuations: 
Spectrum is amplified more at small scales than at large scales. The difference be- 
tween the present study and the experiment may be due to the assumptions made 
in the experimental data analysis, such as turbulence isotropy /homogeneity, and 
negligible pressure fluctuations, which may be too crude in light of the simulation. 
Velocity fluctuation variances are axisymmetric as shown in Sec. 2.1, and thermo- 
dynamic property fluctuations are not isobaric and decay rapidly behind the shock 
wave as shown in Sec. 3.3. The effects of these imperfect assumptions on the data 
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kilk 0 

FIGURE 5. Changes in one dimensional spectra across the shock wave which 
is at k a x, = 12.29: lines for upstream spectra at k 0 X\ = 10.59, and symbols for 

downstream spectra at k Q x\ = 13.68. , • : £ 2 (^ 2 )) 1 * : -^ 1(^2 )i i 

A : £3(^2 )• 

analysis are not clearly documented in Keller et al. (1990). 

Figure 6(A) shows the evolutions of Taylor microscales (A 0 ) and the transverse 
density microscale (A p ), which axe defined as 



respectively. All the microscales decrease significantly across the shock wave: the 
streamwise Taylor microscale by about 50%, the transverse Taylor microscales by 
about 20%, and the density microscale by about 30%. Mach number dependence of 
Taylor microscale change predicted by LIA is shown in Figure 6(B). The higher the 
Mach number, the Taylor microscales are reduced further through the shock wave. 
The reduction is more pronounced in the shock-normal direction. The reduction ob- 
served in the simulation agrees well with the LIA prediction. The Taylor microscale 
which was reported to increase (Debieve et al. 1986) was the time scale, not the 
length scale (Debieve 1992, private communication). If the mean velocity decrease 
across the shock is properly accounted for, their experimental result is consistent 
with the present simulation and analysis. 
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FIGURE 6(a). Evolution of microscales throughout the computational domain 

for Mi = 2.0: Ai, A 2 , A3, A^. Vertical lines denote the 

boundaries of shock intermittency. 



FIGURE 6(b). Change of Taylor microscales and dissipation length scale across 

the shock wave predicted by LIA. Ai, A 2 & A 3 , / £ . Symbols for 

DNS: • A 2 , x / e . 


Integral scale (A/) of turbulent fluctuation /' in the x 2 -direction is defined from 
its two-point correlation, C/(r;x 1), defined as 


C f (r- Xl ) = <)/*(*! '*2 + r ,x 3 ,<) 

/ (^1 1 ^ 3 ^)/ *^3 >0 
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FIGURE 7. Evolution of integral scales throughout the computational domain for 

Mi = 2.0: A Ua , A Ul , A Ua , A p . Vertical lines denote the 

boundaries of shock intermittency. 

where the average is taken over time and homogeneous directions (a?2- and x 3 - 
directions). The integral scale (A/) is, then, defined as 

A/(xi) = / Cf(r;x x )dr, 

Jo 

where the upper limit of the integration is replace by L/2 when dealing with numer- 
ically simulated field with L being the computational box size in the ^-direction, 
where the periodic boundary condition is enforced. Figure 7 shows the evolutions of 
four integral scales throughout the flow field. Three integral scales (A Ul , A Ua , and 
A p ) undergo reductions across the shock wave, most significantly in A U2 by about 
45%, while A« a increases by about 30%. Mach number dependence of the integral 
length scale change can be predicted by LIA. For the shock wave with M\ = 2.0, 
the ratio of the downstream to the upstream integral length scale (with the von 
Karman upstream spectrum) is 0.91, 0.60, 1.46 for A Ul , A U2 , and A Ua , respectively. 
The simulation results agree well with the LIA predictions considering the difference 
in the upstream energy spectrum shape (see Fig. 4). 

Most widely used length scale in turbulence modelling is the dissipation length 
scale (/ € ), defined as 


h = pq 3 /*, 

where e is the dissipation rate of turbulence kinetic energy, which includes contribu- 
tions from both solenoidal and dilatational motions. Figure 8 shows the evolution 
of the length scale l € . The dissipation length scale also decreases across the shock 
wave. Just behind the shock wave, the length scale undergoes rapid increase as 
was the case with the streamwise Taylor microscale (Figure 6(A)), due to the rapid 
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FIGURE 8. Evolution of the dissipation length scale throughout the computational 
domain for M\ = 3.0. Vertical lines denote the boundaries of shock intermit tency. 

decay of the acoustic waves (or, the dilatational motions). The Mach number de- 
pendence of the dissipation length scale change predicted by LIA is presented in 
Figure 6(B). The length scale is reduced for strong shock waves, while it shows a 
mild increase for shock waves with M\ < 1.65. For weak shock waves, TKE and 
its dissipation rate is comparably amplified to give slight increase in l € across the 
shock wave, while TKE amplification saturate much faster than vorticity variance 
amplification to give the reduction in the length scale (Lee et al. 1993). The length 
scale increase observed by Honkan et al. (1992) at M\ = 1.24 (the equivalent shock 
normal Mach number in their experiment is 1.24 not 1.62) might be explained els 
the phenomenon occurring for weak shock waves, but the analyzed experimental 
results are not in quantitative agreement with the simulation and the analysis: LIA 
predicts less than 10% increase, while the analyzed experimental data shows more 
than 600% increase. This difference seems to suggest that the assumptions used 
in the experimental data analysis may be too crude, such as negligence of pressure 
fluctuations and applicability of Taylor’s hypothesis in high intensity turbulence, 
and should be examined carefully. 

2.3 Thermodynamic quantities 

Thermodynamic fields which are obtained from the freely decaying turbulence 
(Lee et al. 1991b) and prescribed at the inflow are nearly isentropic. As the flow 
passes through the shock wave, all the fluctuations are amplified, followed by a 
decay. A general assumption on the relation between thermodynamic fluctuations 
is polytropic (with a polytropic exponent n), where 

p' _ p' _ n T" 
p p n — If 





342 


Sangsan Lee 



FIGURE 9. Evolution of exponents for the fluctuations in thermodynamic variables 

for M\ = 2.0: n pp , ti p t, c p t • Vertical lines denote the boundaries 

of shock intermittency. 


For polytropic fluctuations, specification of one property fluctuation and the poly- 
tropic exponent is enough to describe the thermodynamic fluctuations. Based on 
the above relations, different poly tropic exponents can be defined using normalized 
rms fluctuations ( n pp , n p r) and the correlations between instantaneous fluctuations 
(c pT ) as 


and 



n p T = 1 + 



c p t = 1 + 


p p'T" 
T ~p* 


For weak shock waves with M\ < 1.20, relations between thermodynamic prop- 
erty fluctuations are close to isentropic (n = 7 ) throughout the flow field (Lee et al. 

1993). 

In order to check the polytropy for the strong shock case, the polytropic ex- 
ponents, n pp ,n p Ty defined above are investigated. If the fluctuations are indeed 
polytropic, the two exponents should be the same, which is defined as the poly- 
tropic exponent. The evolutions of the two exponents are shown in Figure 10. The 
exponents are the same upstream of the shock wave with n pp = n p r ~ 7. Down- 
stream of the shock wave, however, they differ significantly with n pp decreasing and 
n p x increasing. Their return to polytropy is very slow. To further investigate the re- 
lation between instantaneous fluctuations, the correlation between the fluctuations 
of density and temperature (p'T") is studied. 
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The evolution of the exponent is also shown in Figure 9. Upstream of the shock 
wave, the exponent is quite close to the 7 (= 1.40). It drops significantly across the 
shock wave, and its further evolution is rather slow. The change in the exponents 
across the shock wave is found to be consistent with the LI A prediction (shown 
in Figure 10). Upstream thermodynamic fluctuations are polytropic (close to isen- 
tropic), and downstream fluctuations are not isentropic due to significant entropy 
fluctuations produced by the shock turbulence interaction. To properly describe the 
thermodynamic fluctuations in strong shock turbulence interaction, specification of 
at least one thermodynamic fluctuation along with two exponents (t. e., n pp and 
ti p t) are required. 


The shock strength effects on thermodynamic fluctuations for wider range of 
shock normal Mach numbers can easily be investigated through the linear analysis. 
In the following, polytropic exponents downstream of the shock in the interaction 
of solenoidal velocity fluctuations with a shock wave is studied. The effects of the 
shock strength on downstream poly tropic exponents are shown in Figure 10. For 
isentropic or acoustic fluctuations, all the exponents are same and equal to the 
specific heat ratio. For entropic or isobaric fluctuations, n pp and c p t become 0 and 
n p r becomes 2. For weak shock waves with M\ < 1.2, thermodynamic fluctuations 
behind the shock can be regarded as isentropic. As the shock becomes stronger 
beyond this limit, the entropy fluctuation behind the shock cannot be neglected, and 
its importance becomes more dominant for the stronger shock waves. The results 
of the poly tropic exponents from DNS are consistent with LIA predictions with 
the values from DNS systematically deviating from the LIA predictions toward the 
isentropic value of 1.4. This may be due to (incompressible) pressure fluctuations 
associated with dilatation- free velocity fluctuations (Sarkar et al 1991), which 
accompany mainly isentropic thermodynamic fluctuations. 


In order to quantify the importance of entropy fluctuations behind the shock wave, 
the contributions of acoustic and entropic fluctuations to the density fluctuation are 
quantified by the linear analysis and also shown in Figure 10. Since the acoustic 
fluctuations and entropic fluctuations are completely decorrelated in the linear limit, 

the relative importance of entropy fluctuations can be expressed as For weak 

shock waves with A/i < 1.2, entropy fluctuations contribute less than 2% to the 
density fluctuations. However, entropy fluctuations become more important than 
acoustic fluctuations beyond Mi = 1.65. 


In summary, thermodynamic fluctuations downstream of the shock wave are 
found to be isentropic for weak shock waves ( M\ < 1.2) and become non-poly tropic 
for strong shock waves, where the importance of entropy fluctuations are com- 
parable to the acoustic fluctuations. The thermodynamic fluctuations cannot be 
modelled using polytropic exponents in this regime. Therefore, modelling effort 
should be made separately for the acoustic fluctuations and entropic fluctuations. 
Zeman (1993) stressed the need for such a separation for the mean thermodynamic 
quantities. 
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FIGURE 10. Downstream poly tropic exponents and entropy fluctuation contribu- 
tion predicted by LI A. n pp , n pT , c pT , Entropy fluctuation 

contribution. Symbols for DNS: • n pp , o ti p t , * c p r- 


3. Future plans 

The varied evolution of thermodynamic variables in Large-eddy Simulation (LES) 
using different formulations (where no explicit removal of aliasing errors is per- 
formed) has not yet been understood. In order to have a reference case where alias- 
ing errors are removed exactly, a specific volume formulation in solving compressible 
Navier-Stokes equations is being pursued. Large eddy simulation of isotropic tur- 
bulence with a shock wave will be performed once the cause for the difference in 
the evolution of thermodynamic quantities is better understood. 

Numerical simulation will be extended for a more practical situation where the 
turbulent boundary layer is subjected to externally imposed strains: a boundary 
layer under rapid expansion and over a compression ramp. 

This work was produced in collaboration with Prof. S. Lele and Prof. P. Moin. 
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